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ABSTRACT 

The instabilities of some spatially and/or time-periodic flows are 
discussed, in particular, flows with curved streamlines which can support 
Taylor-Gortler vortices are described in detail. The simplest flow where this 
type of instability can occur is that due to the torsional oscillations of an 
infinitely long circular cylinder. For more complicated spatially varying 
time-periodic flows, a similar type of instability can occur and is spatially 
localized near the most unstable positions. When nonlinear effects are 
considered it is found that the instability modifies the steady streaming 
boundary layer induced by the oscillatory motion. It is shown that a rapidly 
rotating cylinder in a uniform flow is susceptible to a related type of 
instability; the appropriate stability equations are shown to be identical to 
those which govern the instability of a Boussinesq fluid of Prandtl number 
unity heated time periodically from below. 
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1 • INTRODUCTION 


Our main concern is with the nature of the centrifugal instability of 
time-periodic flows which interact with curved surfaces. It appears at this 
stage that theory and experiment here are in much closer agreement than is the 
case for the instability of flat Stokes layers. In the latter problem it 
seems that the instabilities which exist experimentally occur at such high 
Reynolds numbers that quasi-steady theory applies -and the solution of the full 
time-dependent stability equations, being damped, are irrelevant. The reader 
is referred to the papers by Cowley and Kerczek, which appear in this volume, 
for a detailed discussion of the flat Stokes layer. In contrast to this 
situation, for curved Stokes layers the solutions of the time-dependent 
equation explain a great deal of the available experimental results, whilst 
the quasi-steady solutions of the equations are irrelevant except in the 
inviscid or high wavenumber viscous limits. 

The existence of a Taylor-Gortler instability in the Stokes layers on a 
torsionally oscillating cylinder was demonstrated experimentally and 
theoretically by Seminara and Hall (1976). Since Rayleigh^s criterion does 
not apply to time-dependent flows, there is no simple way of determining 
whether a given time-dependent basic flow is inviscidly unstable. Indeed, 
since Floquet theory can only be used for equations with time-periodic 
coefficients, unless the basic state is varying periodically in time, it is 
not even clear how to define instability. 

The onset of instability predicted theoretically is consistent with the 
available experimental observations; however, at high Taylor numbers it is 
found experimentally that a secondary instability which progressively doubles 
the axial wavelengths occurs. At some stage in this process the flow acquires 
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an azimuthal dependence; there is as yet no adequate theoretical description 
of this regime, and it seems likely that only a fully numerical investigation 
will reveal its structure. 

More recently the instability mechanism found by Seminara and Hall (1976) 
has been shown experimentally and theoretically to occur in a wide range of 
flows of practical importance. These flows vary (periodically) in time and at 
least two spatial directions and exhibit the so-called steady streaming 
phenomenon associated with the work of, for example, Rayleigh (1883), 
Schlichting (1932), Stuart (1966), and Riley (1967). In these flows the 
Reynolds stresses in the boundary layer of the first-order oscillatory viscous 
flow drive a mean motion which may or may not be confined to the Stokes 
boundary layer. In fact, the structure of the steady streaming depends on the 
size of R s , a Reynolds number associated with the mean motion. The 
importance of R g was first explained by Stuart (1966) who showed that for 
R s >> 1 the steady streaming decays to zero in an outer boundary layer of 
relative thickness R+ ^ # it turns out that the instability of the basic 
oscillatory flow also occurs for R^ » 1 and that the steady streaming is 
then driven by both the Reynolds stresses of the oscillatory flow and the 
instability. It is likely that in some situations the instability might cause 
the premature separation of the steady streaming boundary layer. 

The thermal convection analogue of the torsionally oscillatory cylinder 
problem has apparently not yet been investigated. The instability equations 
for an infinite layer of fluid heated sinusoidally from above or below will be 
derived in this paper. It turns out that if the fluid has Prandtl number 
unity then the same equations govern the instability of the flow past a 
rapidly rotating cylinder in a uniform flow. The analogy between the problems 
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is similar to that which is known to exist between the steady rapidly rotating 
Taylor vortex and steady Benard problems. Some preliminary results for the 
solution of the eigenvalue problem governing the instability of these flows 
are given. 

The procedure adopted in the rest of this paper is as follows: In Section 

2 the instability of the flow on a torsionally oscillating cylinder is 

discussed. In section 3 the generalization of this instability to spatially 

varying time-periodic flows is discussed. In Section 4 we discuss the 

analogy between the instability problems for the flow around a rapidly 
rotating cylinder and that driven by the time-periodic heating of a fluid 
layer. Finally, in Section 5 we draw some conclusions. 

2. THE TORSIONALLY OSCILLATING CYLINDER PROBLEM 

Consider the viscous flow induced by an infinitely long circular cylinder 
of radius R rotating about its axis with angular velocity ft cos ait. The 

time-periodic boundary layer on the cylinder is taken to be small compared 
to R so that, 

(-) 1/2 « R, 

and the velocity field u^ is then given by, 

u = fiR(0, V(n, t), 0), 

where 

x = oit, n = {r - R}{^} /2 


and 
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v(n, t) = cos{t - n}e n . 


( 2 . 1 ) 


Thus V(r) , x) is the usual Stokes velocity profile and it is convenient to 
define the Taylor-Gortler number T by 


T = 2 


ft 2 R 

7 71 v V 2 ’ 


( 2 . 2 ) 


and it is assumed from now on that T ~ 0(1). 

The above time-periodic basic state is then perturbed to a disturbance 
which is periodic along the axis of the cylinder with wavelength based on 

2v V k 

(~) 2 the boundary layer thickness. After some manipulation we find that 

the radial and azimuthal velocity components U, V satisfy the coupled 
system: 


£^-2 ~ k 2 }U = 2k 2 TV (p , t)V , 


3n 


(2.2a) 


<£v = /2 U, 

Oh 


(2.2b) 


where 


£ ^ - k 2 - 2 f 

3n 2 5t 


The system (2.2) must be solved subject to 


„-V=fI-0, 1 = o, 


and 


(2.3) 


u = v = o, n = », 


(2.4) 
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so that the no-slip condition is satisfied at the wall and the disturbance 
decays exponentially to zero at infinity. There is, of course, no 
justification for a quasi-steady approximation if T is 0(1) since the time 
dependence of the partial differential system is in no sense slow. 

On the basis of Floquet theory we expect solutions of (2.2) of the form 

U = e aT u(n, t), V = e aT v(n> x) 


where u, v are now periodic functions of x# The Floquet exponent a is a 
function of k and T and must be chosen such that the partial differential 
equations and boundary conditions are satisfied. If we write 

U = l e inT u n (n), V = l e inT v n (n), 

—03 —00 


then the eigenvalue a is determined by the system 


rd 2 ,2 u d 2 .2-, k 2 T r -n(l+i) . -n(l-i) i 

[-T - k - ® ' 2ln ] [-T ' k K * — [ e Vi * e Vil- 


li 

dn" dn 

[^T - k2 - ° - 2in ] v n = 5= t (1 + i)e' n(1+l) u n _l + (1 " i ) e " n(1 " l) VJ’I 


dn 


n 


*(2.5) 


du 


u = v = - — =0, n = o, 
n n dn 


u = v =0, n = °°> 
n n 


for n = 0, db 1 , ± 2 
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The eigenvalue problem (2,5) can be solved numerically by truncating n 
such that |n| <_ N and finding a = cf(N, K, T). Numerical calculations by 
Seminara and Hall (1976) showed that instability occurs for 

T > 164.42 

and the critical value of the wavenumber is k = 0.859. The neutral curve 
calculated by Seminara and Hall is shown in Figure 1. The calculations also 
showed that u n = 0 for n odd and that v n = 0 for n even. These 
results are in excellent agreement with the experiments performed by Seminara 
and Hall (1976) and Barenghi, Park, and Donnelly (1980) who visualized the 
instability using dye visualization. 

However, in both sets of experiments a secondary instability at a Taylor 
number about 30% above the critical was observed. The linear theory of 
Seminara and Hall showed that when instability occurs the function Vq(ti) is 
zero so that there is no mean flow around the cylinder. When the secondary 
instability occurs the vortices initially interact with their neighbours to 
generate larger vortices. Subsequently this larger set of vortices interact 
to produce an even larger set and this process appears to continue without any 
equilibrium state ever being reached. At some stage in this process the 
azimuthal velocity field develops a mean component and there are signs of 
nonaxi symmetry. 

At this time very little is understood theoretically about how and why 
this secondary instability takes place. Seminara and Hall (1977) showed that 
the first mode bifurcates supercritically at the critical Taylor number. Hall 
(1981) showed, using an approximate approach, that an axisymmetric mode with 
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twice the wavelength of the first mode causes the first mode to lose stability 
very close to the experimentally observed critical Taylor number. However, it 
should be stressed that this approximate approach, which is similar to that 
recently used for mode interactions in spatially varying flows, cannot be 
rigorously justified; nevertheless, it would appear that it describes very 
well the first stages in the secondary instability regime. It is of interest 
to note that in the steady Taylor and Gortler vortex instabilities the onset 
of the secondary wavy vortex modes leads to a new equilibrium state. 

The possible role of nonaxisymmetric modes in the secondary instability 
problem is still not understood; certainly it is known from the work of Duck 
and Hall (1981) that 0-dependent modes occur quite close to the axisymmetric 
critical Taylor number and are more stable than the axisymmetric one. 
Moreover, we can show that when nonaxisymmetric modes are present a steady 
streaming phenomenon occurs and the instability can no longer be confined to 
the boundary layer. In this case the velocity field decays to zero in an 
outer layer of relative thickness R^ 2 where R g is a steady streaming 
Reynolds number. 

It would appear that little further theoretical work on the secondary 
instability stage is possible until more precise experimental results become 
available. The flow visualization methods used previously are not sensitive 
enough to determine the detailed flow structure after the secondary 
instability occurs. 
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3. INSTABILITIES IN FLOWS EXHIBITING STEADY STREAMING 

The prototype problem discussed in the previous section shows how a 
centrifugal instability mechanism can occur when a Stokes layer interacts with 
a curved surface. However, in general it is known that when such an 
interaction takes place a secondary steady streaming is set up. The nature of 
the steady streaming depends crucially on a steady streaming Reynolds number 
R s whose importance was explained clearly by Stuart (1966). For large values 
of R s the streaming decays to zero in a steady outer boundary layer of 
relative thickness R g ^ whilst for small values of R g the steady motion is 
confined to the Stokes layer. We shall see below that the instability 
mechanism described in Section 2 occurs for R s » 1 and has an 0(1) effect 
on the steady streaming boundary layer. 

The possibility that the instability described in Section 2 could occur in 
spatially varying oscillatory flows over curved walls was overlooked until 
Honji (1981) investigated the classical steady streaming flow induced by 
oscillating a cylinder transversely along its axis. In addition Honji 
investigated flows over wavy walls and over steps and found that a vortex 
instability of these time-periodic flows occurred at sufficiently large 
amplitudes of oscillation of the fluid velocity at infinity. We will discuss 
only the transversely oscillating cylinder problem in detail; related results 
for wavy walls, pipe flows, and spheres will be mentioned briefly later. 

Ift his experiment Honji (1981) investigated the classical steady streaming 
problem of Schlichting (1932), Stuart (1966), and Riley (1967). For 
sufficiently small amplitude high frequency oscillations, he found that the 
flow remained two-dimensional and was consistent with theoretical 
predictions. However, at a critical frequency dependent amplitude of 
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oscillation the flow became three-dimensional in the neighbourhood of the 
positions on the cylinder where the tangent plane was parallel to the 
direction of motion of the cylinder. At a still higher critical amplitude of 
oscillation, Honji found that the flow became "turbulent and separated." We 
will show below that a linear stability analysis predicts the first critical 
amplitude of oscillation whilst a nonlinear theory suggests that a finite 
amplitude instability might cause the basic steady streaming boundary layer to 
separate prematurely. The details of the calculation outlined below can be 
found in the paper by Hall (1984). 

Suppose that a circular cylinder of radius a oscillates with velocity 
Uq cos rnt along a diameter in a fluid of viscosity v . The parameters which 
govern the two-dimensional flow are 



(3.1a) 


U 

wa 


R 

s 



U)V 


x 2 3< 


(3.1b) 

(3.1c) 


The frequency parameter 0 is taken to be large so that the boundary layer on 
the cylinder is thin compared to its radius whilst the amplitude parameter X 
is taken to be small. Before specifying the size of the steady streaming 
Reynolds number R g we write down the Taylor number T based on (v/w) ^ 
the boundary layer thickness. We obtain 
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T = 


2 3/ 2 „a 


av 


'/2 , 3/2 


■2 3/2 R f 1 ^ 
s 


so that if instability occurs at 0(1) values of T we must take R g 

formally to be 0 ( g ^ ) • thus we take the limit 0 » with R g g~ ^2 held 

fixed. We further note that in this limit X ~ g ^ so that the boundary 
layer on the cylinder is essentially a Stokes layer. In Figure 2 we have 
sketched the cylinder and the various regions of interest. The inviscid slip 
velocity for the basic flow has maxima at 9 = ± tt/ 2, and since the curvature 
of the boundary is constant we expect that these will be the least stable 
locations. A WKB approach to the instability problem quickly shows that 
0 = ± ir/2 are turning points of the expansion. In fact, they are second- 
order turning points so that an inner region of angle 0(g~*^) is needed at 
0 = ± -it / 2 . 

Suppose then that the basic flow is perturbed to a disturbance periodic 
along the axis of the cylinder with wavelength The linear stability 

1C 

equations can be reduced to 


2 p 2"^ 4 sin 0 v T ^ 

L(— ■ — ■j - k^)U = 2k 2 T sin 8 v Q V - |£ + 0(g" ^ ) , 


9n 


0 




(3.2) 


LV = 4 sin 0 — 2 . U + 0(0 ^ ), 
an 


L = 




2 5/4 T 1/2 sin 0 v 




0 9_ 
90 


where 
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wit h 


Vq = cos t - cos(t - ri)e r '. 


(3.3) 


The Taylor number T is expanded as 


T 




+ ... 


whilst near 8 = ir/2 the velocity components U, and V expand as 


1 / 8 

u = u Q (n, t , $0 + 8 u^n, t, $) + ••• 
v = v 0 (n, t, $) + 8 1 ^ 8 v 1 (n, t, $) + ••• 

1 /8 

where $ = 0 (0 - n/2). The functions Uq, Vq can be written 


(Uq, V q ) = A($)(u 0 (n, x), v Q (n, t)) 


where A is an amplitude function to be determined whilst (uq, vq) satisfy 



21(2 T o v o V 



3v r 


ffK = 4 v 


3u^ 

u 0 = Tn" 


= v 0 = °, 


(3.4) 


u o> v o + o* n ”* 


n = 0, 
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This eigenvalue problem was solved by Hall (1984) in exactly the same manner 
as that described in Section 2. The neutral curve k = k(Tg) is shown in 
Figure 3; we see that instability is predicted locally near 9 = tt/ 2 for 


T 0 > 11.99. 


(3.5) 


The amplitude function A($) is found at higher order to satisfy 


2 

~2~ + y[T, - $ 2 ]A = 0. 
d$ 


(3.6) 


Here y is a positive constant which must be calculated numerically; the 
solutions of this equation, which decay when $ ->■ ± 77/2 are 


A = U n (-n - 1/2 , 2y 1/4 $) 


where U n is the nth parabolic cylinder function and 


T = T = - f n ]/Z 1 

1 r v - 


The most unstable mode corresponds to n = 0 and gives 


A = e " u ^ $2/2 
A 0 > 

with 


ho -~p - 5 - 51 - 

VU 


Thus the critical value of R g is 
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R 

s 


= R 


sc 


4.24[B^ 2 + .460 ^ + 


■]. 


(3.7) 


or in terras of X the critical oscillation amplitude corresponds to 


X = X 

c 


2.06 

Y4 


[i 


.23 


A 




(3.8) 


In Figure 4 we have plotted the results predicted by (3.8) and Honji's 
results; there seems little doubt that the instability mechanism is 
responsible for the onset of the three-dimensionality described by Honji. 

The development of this instability into a weakly nonlinear region leads 
to some surprising results. We choose to write the Taylor number T in the 
form 

T = T a + 0~ ^ T 
Uc 


where Tq c is the critical value of Tq. A Stuart-Watson expansion with the 
fundamental terms of order shows that the generalization of (3.6) into 

the nonlinear stage is 

2 

+ y [f - $ 2 ]A = fA 3 . (3.9) 

d $ 2 

Numerical calculations by Hall (1984) showed that y is negative so that 
finite amplitude solutions bifurcate subcritically at the critical Taylor 
number. The determination of y requires the solution of several partial 
differential systems so the possibility that a numerical error was made should 
not be overlooked. However, if y is indeed negative then it suggests that 
experimentally three-dimensionality could be induced by sufficiently large 
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subcritical perturbations. The nonlinear theory shows that the steady 
streaming boundary layer is, at least in an 0(g _1/8 ) neighbourhood of 
6 = it/2, driven by both the instability and the basic oscillatory flow. If we 
define an outer boundary layer variable £ by 


5 - i~) 


,6 *A 


then the outer steady boundary layer is given by 


with 


— v 3 / 8 — v Va 

v V u = — 8 




which must be solved subject to 


(3.10) 


<l> = 0, 5 = 0 


'P, * 0, 5 + 


(3.11) 


• 23 A < / 2 a^], { - 0 


If we set A - 0 above we obtain the equations governing the attachment of a 
steady streaming boundary layer within a g neighbourhood of 0 = n/2. 

In that case (3.10) is solved subject to 


^ = 0 , $ = 0 , 


(3.12) 
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and the symmetry of (3,9) about $ = 0 means that (3*12) can still be applied 


since either 


dA 

A or - = 0 at $ = 0. We note that for large $ the 


, £ + 0> so that, assuming that the 


3 $ 

condition (3.11) reduces to \l> + T~ 0 

£ uc 2 ^'^ 

boundary layer remains attached for finite values of $, the extra term 

dA 

proportional to A — in (3.11) merely produces an origin shift in the large 


$ asymptotic solution of (3.10). The linear eigenfunctions A r ($) for 
n t 0 all have intervals where A n ($)A^($) is positive so that the 
possibility exists that in the nonlinear regime the slip velocity in (3.11) 
might change sign at some value of $• If the magnitude of the inviscid slip 
velocity is sufficiently large where this occurs, then the attached flow 
strategy fails and the steady streaming boundary layer will detach prematurely 
from the cylinder. This possibility does not occur for more general flows 
where the point of attachment of the steady streaming layer and the most 
unstable position do not coincide. In this case the steady streaming driven 
by the instability is weak compared to that of the basic flow. 

For more general steady streaming flows Papageorgiu (1985) has shown that 
both concave and convex curvature lead to this local vortex type of 
instability. Consider then an oscillatory viscous flow adjacent to the wall 
y = 0 induced by a outer potential flow with slip velocity U(x)cos wt. 
If the wall has radius of curvature aR(x), we write down the local Taylor 
number T A defined by 

a J 




9 3/2 2 „2, 

2 U Q U (x) 


av 


Vi Jn 


0 ) 


R(x) 


The mechanism described above for the circular cylinder problem occurs near 

x = x if 
m 
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t ( x„> ~ 11 •" 

T t-°' x- 


Hence with variable curvature the most unstable locations do not necessarily 
occur where the inviscid slip velocity has a maximum; this produces some 
interesting results. Thus for an ellipse oscillating transversly there can be 
two or six local regions of instability depending on the angle of attack and 
the eccentricity of the ellipse. 

Papageorgiu has shown that for concave curvature the corresponding 
condition is 

|Tj = 7.104, 


so that a locally concave wall is more susceptible to instability than is a 
convex one. The result that either a concave or a convex wall lead to 
instability is quite different than the classical results for steady flows; 
however, for a time-dependent flow Rayleigh's criterion does not apply so 
there is no reason why this should not be the case. 

The implications of Papageorgiu's calculations for flows in curved pipes 
and over wavy walls are important. For a curved pipe he has shown that an 
oscillating azimuthal pressure gradient leads to instability at both the inner 


and outer bends; however, the outer bend is the most unstable location. This 
instability leads to rapid local variations in the shear stress at the wall 
which might have important consequences for aortic blood flow. It should also 
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be pointed out that in this problem the turning point structure changes from 


that found 

for 

the 

circular cylinder problem; 

the 

technical 

problems 

associated 

with 

this 

change were overcome using 

the 

expansion 

procedure 


devised by Soward and Jones (1981). A similar difficulty arises when the 
instability of the flow due to a torsionally oscillating sphere is considered; 
in this problem the instability is localized near the equation. As yet the 
generalization of the expansion procedure of Soward and Jones into the 
nonlinear regime has not been given. 


4. THE LINEAR STABILITY EQUATIONS FOR THE FLOW AROUND A RAPIDLY ROTATING 
CYLINDER 


Consider the flow of a viscous fluid of kinematic viscosity v around a 

circular cylinder of radius a. The cylinder rotates with angular velocity 
V 

— whilst the flow a long way from the cylinder has sped U in the x 
direction. We follow the notation of Moore (1957) and define the 
parameter e by 


e 



(4.1) 


and a Reynolds number R by 



(4.2) 


Moore has discussed the above flow in the limit' e -*• 0; here we will also 

assume that R » 1. For large values of R it is known from the work of 

Glauert (1957), Moore (1957), and Wood (1957) that a boundary layer of 
-l/ 9 

thickness aR z is set up on the surface of the cylinder. Here we will 
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consider the possible instability of this boundary layer to a Taylor-Gortler 

vortex type of perturbation. An examination of the effective Taylor-Gortler 

number for such a flow suggests that instability is likely for 0(1) values 
of the parameter 


T = 



(4.3) 


hence we limit our investigation to the limit e * 0 with T, which we will 
refer to as the Taylor number, held fixed. On the assumption that e ~ R~ 
we can show from the work of Moore (1957) that in the boundary layer u and 
v, the radial and azimuthal velocity components, are given by 

u = VeR~ fu(n, 0) + 0(R" 1/ 2)} > ( 4#4a) 


v - V{1 + ev(n, 0) + 0(R -1 ) } . 


(4.4b) 


Here n is a boundary layer variable defined by 


-- H Ji) (!) 1/2 


(4.5) 


whilst u and v can be obtained from Moore's work; in fact, we shall only 
need an explicit form for v which is given by 


v 


£e i0 ~h( 1+ i) 


+ complex conjugate. 


(4.6) 


Thus v is just the Stokes layer velocity profile with 0 taking the place 
of the time variable. It is this correction to the irrotational flow v = l/ r 
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which leads to instability. This result follows from the fact that the 

irrotational flow v ~ 1/r is neutrally stable according to Rayleigh's 

criterion. 

Suppose that we define the dimensionless variables t and Z by 




where z represents distance along the axis of the cylinder and t denotes 
time. The flow given by (4.4) is perturbed by writing 

R~ 1/2 [u + 0(R _1/2 )] + U(n, 0, T)e lkZ , (4.7a) 

^ - [1 + v(n, 0) - n + 0(R h] + V(n, 0, r)e lkZ , (4.7b) 

^ = W(n, 0, z)e ikZ . (4.7c) 


Thus k is the wavenumber of the vortex type of instability and, unlike the 
usual case with a centrifugal instability mechanism, the three components of 
the distrubance velocity field are of comparable size. The disturbed velocity 
field (4.7) can be substituted into the Navier-Stokes equations, and after 
some manipulation we obtain the linear stability equations: 

3n 

(72 - k2 - 2 le - 2 !r> v - ' 2 ™ !!• 

3r| 


(4.8b) 
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• 2 !e - 2 It!" - lkp «^- 


(4.8c) 


Here P is the pressure perturbation scaled on pvR~ ^ 2 v 2 and terms of 
relative order R 2 have been neglected. The continuity equation 
corresponding to (4.8) is 


au 

9n 


+ 


ikW 


0. 


(4.9) 


It is convenient at this stage to eliminate P and W from (4.8) and (4.9) 
to give the coupled pair of equations: 




2 " k l u - ^ 2 V, 

9n 



(4.10) 


which must be solved subject to 



= 0, 


n = o, oo. 


Thus we assume that the instability is confined to the boundary layer; this is 
to be expected so long as we ignore nonlinear effects. 

Suppose next that a viscous fluid of Prandtl number unity and viscosity v 
occupies the region z 0 and that the plane z = 0 has temperature 
2 Tq cos mt. We define dimensionless variables t and q by 


x = rnt, 


(4.11a) 
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n 



(4.11b) 


so that the motionless state has temperature field T given by 


T = T q v(t, n) 


(4.12) 


where v is as defined by (4.6) with 0 replaced by x. Thus the 
temperature is zero except for a thin region of thickness A = ^ whilst 
in the boundary layer the temperature is given by the Stokes layer velocity 
profile. In the usual way the basic state is now perturbed by a disturbance 
with horizontal wavenumber k scaled on A If the Boussinesq 
approximation is made then, following Chandrasehkar (1963), Chapter II, we 
obtain 

- k 2 }{^~Y - k 2 - 2 -|— }W = RQ, (4.13a) 

3n 8n 



(4.13b) 


where 0 is the temperature perturbation and W the z velocity component 
of the perturbation. The parameter R appearing in (4.13) is the Rayleigh 
number defined by 

R = gaA^ Tq v ^ (4.14) 


where a is the coefficient of volume expansion. If the plane n = 0 is 
rigid then (4,13) are to be solved subject to 


9 . „ . ™ , 0 , * - 0 , 


(4.15) 
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Thus if we seek neutral solutions of (4.10) with ~ = 0 we see that the 

oT 

resulting eigenvalue problem is identical to that given by (4.13), (4.15) if 
we associate T with R/4/2 and 0 with t respectively. Thus there 
exists an analogy between the neutral stability problems for the flow past a 
rapidly rotating cylinder and a time-periodically heated fluid layer. The 
analogy corresponds to that between the steady Taylor vortex instability 
between circular cylinders rotating at almost the same speed and Benard 
convection between rigid walls. The eigenvalue problem (4. 13)— (4.14) can be 
solved using the method outlined in Section 2; so far we have found one mode 
of instability with a minimum value of R ~ 54. However, other modes exist, 
but it is not yet clear which is the most unstable. 


5. CONCLUSION 

The instability mechanism found several years ago for the flow around a 
torsionally oscillating cylinder exists in many flows of practical 
importance. When the basic oscillatory flow which supports the instability is 
spatially varying, the instability becomes* spatially concentrated near the 
most unstable positions. Moreover the instability drives the secondary steady 
streaming boundary layer at the same order as does the basic oscillatory 
flow. Recently it has been shown by Hall (1985) that this mechanism can lead 
to the instability of high frequency Tollmien-Schlichting waves interacting 
with a wall of either convex or concave curvature. 

There are many flows where the basic state has a steady component of the 
same order as the oscillatory part, in such flows there exists the possibility 
of interactions between the instability mechanisms associated with the steady 
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and unsteady components. An obvious example of this is the Taylor vortex 
problem with a fixed outer cylinder and an inner cylinder rotating with 
angular velocity Q[ 1 + ecos wt]. For e « 1, the flow is susceptible to 

steady Taylor vortices at sufficiently large vaules of fi whilst for e » 1, 
and a) » 1 the Stokes layer mode is possible. For intermediate values 
of e the modes interact in an as yet undetermined manner; the ratio of the 
axial lengthscales for the different modes is large so that some progress 
could probably be made asymptotically. 

Another interaction problem yet to be investigated is that between the 
Tollraien-Schlichting and Taylor-Gortler modes of instability. Experimentally 
it is known that local transition in the flat Stokes layer takes place as 
Reynolds numbers as low as 300. It follows from (2.2) that for Stokes layers 
on curved walls with local radius of curvature R such that 



the centrifugal mode will occur first. This effectively means that in most 
situations the mode is likely to be dominant. 

The preliminary results, which we have obtained so far for the 
sinsuiodally heated fluid layer, suggest that a convection mode of instability 
similar to that found by Serainara and Hall (1976) exists for that flow. As 
yet no experimental investigations of this problem have been made; it will be 
interesting to see how this flow evolves at high Rayleigh numbers. 
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